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Abstract: Wind erosion is one of the main drivers of soil loss in the world, which affects 20 million 
hectare land of Iran. Besides the soil loss, wind erosion contributes to carbon dioxide emission from the 
soil into the atmosphere. The objective of this study is to evaluate monthly and seasonal changes in 
carbon dioxide emission in four classes i.e., low, moderate, severe and very severe soil erosion and the 
interactions between air temperature and wind erosion in relation to carbon dioxide emission in the 
Bordekhun region, Boushehr Province, southwestern Iran. Wind erosion intensities were evaluated using 
IRIFR (Iran Research Institute of Forests and Ranges) model, in which four classes of soil erosion were 
identified. Afterward, we measured carbon dioxide emission on a monthly basis and for a period of one 
year using alkali traps in each class of soil erosion. Data on emission levels and erosion classes were 
analyzed as a factorial experiment in a completely randomized design with twelve replications in each 
treatment. The highest rate of emission occurred in July (4.490 g CO2/(md)) in severely eroded lands 
and the least in January (0.086 g CO2/(m*d)) in low eroded lands. Therefore, it is resulted that increasing 
erosion intensity causes an increase in soil carbon dioxide emission rate at severe erosion intensity. 
Moreover, the maximum amount of carbon dioxide emission happened in summer and the minimum in 
winter. Soil carbon dioxide emission was just related to air temperature without any relationship with soil 
moisture content; since changes of soil moisture in the wet and dry seasons were not high enough to 
affect soil microorganisms and respiration in dry areas. In general, there are complex and multiple 
relationships between various factors associated with soil erosion and carbon dioxide emission. Global 
warming causes events that lead to more erosion, which in turn increases greenhouse gas emission, and 
rising greenhouse gases will cause more global warming. The result of this study demonstrated the 
synetgistic effect of wind erosion and global climate warming towards carbon dioxide emission into the 
atmosphere. 
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1 Introduction 


Accelerated erosion due to human activities, can detach the light and fine particles of soil. Wind 
erosion occurs when forces are exerted by wind, and overcoming the force is necessary for 
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holding particles on the surface of the ground. Wind erosion is the most common form of soil loss 
in arid and semi-arid areas. Wind and water erosion leads to the elimination of surface vegetation 
cover, which has a great effect on soil carbon storage or release (Montgomery, 2007; van Oost et 
al., 2007; Quinton et al., 2010). More than two-thirds of lands in Iran are arid and semi-arid, which 
are susceptible to different types of erosions. Wind erosion is a common phenomenon in arid and 
semi-arid areas of Iran (Modarres, 2008; Amiraslani and Dragovich, 2011). 

During recent years, the intensity of erosion has increased and leads to dust storms, and has 
negative effects on industry, economy, urban and rural life. Under global warming climate 
condition, the precipitation pattern changed and the intensity of wind erosion increased. It is 
estimated that soils contain more than 80% of biosphere carbon in the Earth (2376-2450 Pg C in 
the 0-2 m soil depth) (Chappell et al., 2016). Unreasonable human activities such as abused land 
use, intense livestock grazing, abandoned lands and destructed vegetation cover have caused a 
great acceleration in soil erosion (Giorgi, 2006; Chappell et al., 2016). Soil erosion is one of the 
main sources of carbon dioxide emission into the atmosphere (Lal, 2019). 

As reported, 549x109 hm? of the Earth's surface is affected by wind erosion, of which 296x109 
hm? is severely affected (Lal, 2003). Wind erosion is known as one of the greatest threats to 
sustainable socio-economic development and it is one of the greatest causes of land degradation in 
arid and semi-arid areas of the world. It is estimated that 46.696 of the Earth degraded by wind 
erosion (Peijun et al, 2004). There is a direct relationship between soil degradation and 
desertification with loss of biodiversity, threat to food security, increase in poverty and risks of 
accelerated greenhouse effect and climate change (Bennett, 1939; Jacks and Whyte, 1939; UNEP, 
1992; Oldeman, 1994; Lal, 1998; Scherr, 1999). 

The destiny of carbon dislocated during the erosion process is the subject of scientific 
discussions. There are two main views: (1) soil erosion acts as a carbon sink supported by 
sedimentologists. They claim that large amounts of carbon may be sequestrated through buried 
organic matter and sediments transmitted by various erosion processes. They hint that oxidation 
rates are lower in sediments than in original soils, due to the deep burial, reaggregation and 
physical protection of dispersed materials (Renwick et. al., 2004; van Oost et al., 2004). And (2) 
soil erosion acts as a carbon source, and soil scientists and agronomists support this idea. They 
argue that breaking down soil aggregates causes great emission of carbon into the atmosphere; 
also the rate of mineralization is higher in sediments than in original soils, which causes the 
depletion of soil organic carbon pool (Lal et al., 2004). Due to the significant links among soil 
erosion, carbon content and anthropogenic activities during land use change, it is beneficial to 
accurately evaluate the carbon transfers associated with soil erosion (Kuhn et al., 2009). 

Temperature is another important factor influencing carbon emission into the atmosphere (Lu et 
al., 2013; Zhang et al., 2016; Zhao et al., 2017). Increasing atmosphere temperature leads to 
terrestrial ecosystems response in the case of greenhouse gasses fluxes (Oechel et al., 2000; 
Pereira et al., 2013). For example, Zhu et al. (2015) stated that climate warming increased 
seasonal emission of carbon dioxide and uptake of CH4. However, Dijkstra et al. (2013) found that 
climate warming did not affect GHGs (greenhouse gases) emission or uptake. These differences 
are generally related to the variations in soil properties such as texture, moisture and nutrient 
content (Johnson et al., 1998; Lagomarsino et al., 2016). 

The organic matter content of the soil is low in arid areas, but low density of SOC and its 
concentration near soil surface make it vulnerable to wind erosion. Most of carbon transported 
with sediment is the labile fraction that is easily mineralized. Therefore, wind erosion-related 
emission of carbon is probably high in these areas. In this study, we tested the following 
hypotheses that stated: (1) different severities of wind erosion might have different effects on soil 
carbon dioxide emission; (2) there might be seasonal variations in soil carbon dioxide emission 
from different landscapes; and (3) there might be a positive relationship between air temperature 
and soil carbon dioxide emission. 


2 Materials and methods 


2.1 Study area 


The study area located in the Bordekhun region, Dayer County, Boushehr Province, Iran 
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(27°59'-28°13'N, 51?28'-51?37'BE; Fig. 1), with an area of about 21.411 km?. The annual mean 
temperature and mean annual precipitation in 20 years are 31?C and 155 mm, respectively. Over 
6496 of precipitation is concentrated in autumn and winter. The eastern wind accounts for 7096 of 
daily observations and the mean speed is 4.73 m/s. Weather data for the study area were obtained 
from the meteorological station of Dayer County, southwestern Iran. 
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Fig. 1 Sampling points of carbon dioxide emission in the Bordekhun region, Iran 
2.2 Determination of wind erosion severity 


This research was conducted in two phases: (1) determination of wind erosion severity using 
IRIFR, an experimental model developed by Ekhtesasi and Ahmadi (1998), and (2) measurement 
of carbon dioxide emission influenced by soil erosion. The main reason for the selection of IRIFR 
is that it is the basis for the executive organization's decisions for wind erosion control in Iran. By 
using IRIFR, it is possible to determine nine geomorphological factors involving in wind erosion 
(Tables 1 and 2). The nine geomorphological factors were scored and rated in this model. Scoring 
each factor was done in each geomorphological facies (Table 1). To determine the severity of wind 
erosion, we calculated the algebraic sum of the scores obtained in each geomorphological facies. 
Table 2 shows the range of scores that can finally determine the displacement of soil particles 
caused by wind in each facies. Based on the regional map of geomorphology, we estimated the 
circumference of each facies specified and the effect of nine erosive factors. Proportionating 


Table 1 Definition of class and estimation of sediment discharge of wind erosion using IRIFR (Iran Research 
Institute of Forests and Ranges) model 


Sediment discharge 


Wind erosion class Rate of erosion Sum of calculated score (thm?) 
I Slight «25 «2.5 
I Low 25-50 2.5-5.0 
III Medium 50-75 5.0-15.0 
IV Severe 75-100 15.0-60.0 


M Very severe 7100 260.0 
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Factor 


Petrology (P) 


Land morphology 
(LM) 


Velocity and 
quality of the wind 
(w) 


Soil and its surface 


(S) 


Density of 
vegetation (V) 


Erosive effect of 
soil surface (E) 


Soil dampness (D) 


Type and 
dispersion of 
aeolian deposits 
(AD) 


Land use 
management (LU) 


Table 2 Method of scoring factors of wind erosion using IRIFR model 


Indicator Score range 
Hard igneous rocks; quartzite; massive limestone and granite 0-2 
Granular texture rocks; hard limestone rock and sandstone 2-4 
Marl and clay; medium-grained alluvium; shale and 47 
conglomerate 
Fine-grained alluvium; coastal sands; aeolian deposits and 7-10 
clay plains 
Mountain 0-2 
Hill 2-4 
Pediment 4-7 
Flat plain 7-10 
Average wind velocity at 10 m, less than 4.5 m/s and no dusty 0-5 
severe winds Z 
Average wind velocity at 10 m, between 4.5 and 5.0 m/s and 5-10 


no dusty sever winds, but can make dust 


Average wind velocity at 10 m, between 5.0 and 5.5 m/s; 
occurrence of one dust storm in the year and 10-15 
1—5 d of sand storm in the year 


Average wind velocity at 10 m, more than 5.5 m/s and dusty 


severe winds Ea 
More than 60% of the soil surface is covered by coarse pebble, 5-09 
cemented or crusted 
About 40%-60% of the soil surface is covered by coarse to 0-5 
medium pebbles 
Less than 4096 of the soil surface is covered by fine pebbles 5-10 
and topsoil texture is sandy-loam or loam 
No pebble is found on the soil surface and topsoil texture is loam 10-15 
to sandy 
Vegetation cover more than 50% and uniform plant 5-0 
distribution 
Vegetation cover between 30% and 50% 0-5 
Vegetation cover between 10% and 30% 5-10 
Vegetation cover less than 1096 10-15 
No sign of wind erosion is found on the surface soil 0-5 
Limited signs of wind erosion are found on the surface soil 5-10 
Relatively extensive signs of wind erosion are found in the 10-15 
surface soil 
Extensive signs of wind erosion are found on the surface soil 15-20 
Soil is always wet under the influence of groundwater -5-0 
Surface soil is wet at some times of a year 0-4 
Surface soil is wet occasionally, but dries rapidly 4—7 
Surface soil is completely dry with rapid drainage 7-10 
Signs of wind deposits are not seen as sand dunes 0-2 
Signs of wind deposits are seen as semi-active nebkha dunes 2-4 
Wind deposits like active and non-active sand dunes or ripple 4-7 
marks are seen in the area 
Active sand dunes, nebkhas and ripple marks are seen in the 7-10 
area 
Sparse forest and rangelands with over exploitation and 
agricultural lands with less than 3 months fallow without a 0-5 
windbreaker 
Rangelands with intense overgrazing 5-10 


Bare lands and deserts having sparse or no vegetation cover 10-15 


Overall score 


—5-15 


—5-15 


—5-10 
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to the amount of force of each factor and its effect on sediment yield, we dedicated a score (Tables 
1 and 2) to it. The sum of the nine scores determines the severity of wind erosion. Factors include 
petrology, land morphology, and their topography, velocity and quality of the wind, soil 
characteristics and topsoil conditions, density of vegetation, erosive effects of soil surface, soil 
dampness, type and dispersion of aeolian deposits and the land use management. After 
determining the wind erosion classes using IRIFR method, 12 sample points were randomly 
selected in each class, in which carbon dioxide emission was measured. 


2.3 Measurement of carbon dioxide emission 


Carbon dioxide measuring instrument was made from polyvinyl chloride (PVC) columns, using 
airtight silicone seal as the lid (Fig. 2). Columns were inserted into the soil in 5 cm depth (Suman 
et al., 2009) in each treatment. Measurement of carbon dioxide emitted from the soil was 
conducted in each treatment using these columns for a year. Carbon dioxide emission sampling 
was initiated in mid-March 2015 and continued monthly till mid-February 2016. A plastic cup 
with 10 mL of 1 N sodium hydroxide (Alkali trap), was placed for one day in each column and the 
top of the column was covered with paraffin film (airtight seal). The trapped carbon dioxide was 
determined by titration of the residual NaOH with 0.1 N HCl after the addition of excess BaClz 
and few drops of phenolphthalein indicator. The same procedure was done for the control 
treatment, using a capped column. Amount of carbon dioxide (mg) collected in the alkali traps 
was calculated using Stotzky (1965) equation: 

Amount of CO2=(B-V)xNxE, (1) 
where B and V are the volumes (mL) of acid used for titrating NaOH of control and treatment in 
the alkali traps, respectively; N is the normality of the acid used; and E is the equivalent weight (g), 
which is 22 g for calculating carbon dioxide content. 


Fig.2 Photo for collecting carbon dioxide emission in the Bordekhun region, Iran 


2.4 Measurements of soil moisture and temperature 


Soil samples were simultaneously taken from adjacent points of carbon dioxide emission 
sampling. Gravimetric soil water content was monthly recorded, and air temperature data were 
obtained from the meteorological station of Dayer County, southwestern Iran. 


2.5 Data analyses 


The factorial experiment in a completely randomized design with twelve replications in each 
treatment (erosion intensity) using SPSS 17.0 was performed to analyze the association between 
soil carbon dioxide emission levels and erosion classes. Analysis of variance (ANOVA) in the 
generalized linear modelling (GLM) procedure was employed to test the differences in the carbon 
dioxide emission in response to air temperature in each wind erosion class. Differences between 
means of carbon dioxide emitted and wind erosion classes confirmed by Duncan test (P<0.05). 
Due to the homogeneity and normality of the data, the Pearson's correlation coefficient was used 
to measure the strength of association between temperature or soil moisture and carbon dioxide 
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emission in different wind erosion intensities. 


3 Results 


We developed Table 3 according to the scoring and rating of the nine factors affecting wind 
erosion. Within this table, the erosion classes of the 16 geomorphological facies in the region were 
defined. Classification of wind erosion intensities in the study area is presented in Figure 3, which 
shows four classes of wind erosion severity. As shown in Table 1, the model can divide the region 
into five erosion classes. In the study area, we couldn't find erosion class I, the minimum erosion 
severity, because of the harsh condition of the study area. Generally, the larger wind erosion 
classes located in areas covered by water erosion deposits, and the minor erosion classes found 
where rocks and hard beds are presented. 


Table3 Evaluation of wind erosion using IRIFR model 


No. Geomorphologic facies P LM W S V E D AD LU SS Eua 
1 Rock outcrop with debris 4 4 4 -2 8 3 4 3 -2 26 II 
2 Coarse debris 6 4 8 -4 15 3 4 2 -3 35 II 
3 Water erosion 5 4 8 0 8 2 6 4 -2 35 II 
4 Irregular slope 6 5 14 5 0 4 6 5 7 52 III 
5 Gully erosion 9 9 18 14 7 18 5 7 14 101 V 
6 Sand beds 10 10 20 15 15 10 10 10 5 105 V 
7 Alluvium of fine grain 9 7 20 4 14 8 8 7 5 82 IV 
8 Sand dunes 10 10 20 15 15 10 9 10 5 104 V 
9 Argillaceous area 10 10 20 10 12 10 10 10 0 92 IV 

10 Side erosion 8 6 15 3 8 6 5 4 0 55 III 

11 Rock outcrop 7 4 10 4 7 4 6 4 -2 44 II 

12 Alluvial fan 8 8 19 10 3 8 9 5 5 75 Ill 

13 Wind erosion 7 10 19 10 11 10 5 5 6 83 IV 

14 Alluvium of coarse grains 9 10 18 15 14 9 9 7 10 101 V 

15 Slight rill erosion 10 9 18 15 13 10 10 7 11 103 V 

16 Alluvium of medium 7 9 17 9 10 7 10 5 5 79 IV 


grains 


Note: P, petrology; LM, land morphology; W, velocity and quality of the wind; S, soil and its surface; V, density of vegetation; E, 
erosive effects of soil surface; D, soil dampness; AD, type and dispersion of aeolian deposits; LU, land use management; SS, sum of 
scores; Eqiass, erosion class. 


3.1 Effects of erosion intensities on carbon dioxide emissions in different months 


The results showed that the amount of carbon dioxide emitted in different months and erosion 
intensities are significantly different. Rising air temperature significantly increased the carbon 
dioxide emitted from the soil in each class of erosion intensity (Table 4). Soil carbon dioxide 
emission was significantly affected in warmer months by increasing erosion intensity. The highest 
amount of carbon dioxide emission occurred in severe erosion intensity areas in July (4.490 g 
CO»/(m?.d)), and the least in low erosion intensity in January (0.086 g CO»/(m?.d)) (Fig. 4). 
Having the same trend, seasonal changes had a significant effect on the carbon dioxide emitted 
from the lands with different erosion intensities. The maximum amount of soil carbon dioxide 
emitted in the summer and the minimum amount in the winter. 


3.2 Correlations between soil moisture and air temperature with carbon dioxide emission 


A positive correlation was observed between air temperature and carbon dioxide emission in all 
erosion intensities. However, there was no significant relationship between soil moisture and 
carbon dioxide emission in any erosion classes (Table 5). Annual precipitation at the study year 
was 185 mm and mean temperature was 27.9°C. The maximum and minimum temperatures 
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Fig. 3 Wind erosion classes using IRIFR model. II, low erosion; III, medium erosion; IV, severe erosion; V, very 
severe erosion. 


Table 4 Result of ANOVA for carbon dioxide emission in different months and wind erosion classes 
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Source Sum of square df Mean square F Sig. 
Corrected model 346.228 47 7.367 185.497 0.000 
Intercept 375.870 1 375.870 9464.806 0.000 
Month 246.564 11 22.415 564.431 0.000 
Erosion 43.937 3 14.646 368.790 0.000 
MonthxErosion 55.727 33 1.689 42.523 0.000 
Error 20.968 528 0.040 
Total 743.066 576 
Corrected total 367.196 575 
6 WLow Bi Medium B Severe li Very severe 


Soil carbon dioxide emission 
(g CO /(m*d)) 


Mar Apr May Jun Jul Aug Sep Oct Nov Dec Jan Feb 
Month 


Fig.4 Soil carbon dioxide emission in different months and wind erosion classes (low, medium, severe and very 
severe) 
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occurred in July (34.9?C) and December (19.2?C), respectively. The maximum precipitation was 
found in December (118 mm) and no precipitation was taken place from April to September. This 
trend in temperature and humidity, in addition to the light texture of the soil, results in surface soil 
dryness. So that the highest moisture content of the soil in December was 1796, and in other 
months varied from 396 to 1396. 


Table 5 Pearson's correlation coefficients of wind erosion classes with air temperature and soil moisture 


Wind erosion class Air temperature Soil moisture 
Low 0.686* —0.534 
Medium 0.629* —0.545 
Severe 0.793" —0.585 
Very severe 0.711" —0.503 


Note: * and ™ indicate significances among different wind erosion classes at P<0.05 and P<0.01 levels, respectively. 


4 Discussion 


Two distinct results were achieved in line with the objectives of this study. The first is that, with 
increasing erosion intensity, the carbon dioxide emission rate increased. This increase in carbon 
dioxide emission was only recorded in areas with severe and very severe erosions. At low erosion 
intensities, carbon dioxide emission was not significantly different from the control. 

Various studies at different regions of the world showed that erosion has a great effect on soil 
carbon loss (Leifeld et al., 2011; Ran et al., 2013), which is in accordance with the current study, 
highlighting the increase of carbon emission according to soil erosion accretion. Soil erosion 
moves topsoil and causes the soil carbon released into the atmosphere. High amounts of organic 
carbon stored at the surface layers of soils, implies the importance of this erodible layers 
conservation against the factors that cause the loss of nutrients directly (i.e., erosion processes) or 
indirectly (via providing the condition of decomposition of soil organic matter). So using some 
methods to reduce wind erosion is beneficial for conserving carbon besides other advantages. 

One of the major causes of erosion in the world is land use change, which influences the 
vegetation and soil of an area and afterward changes the amount of carbon captured in land 
(Houghton and Goodale, 2004). Land use change in Iran is very severe and, among other factors, 
causes the erosion rate of 7 to 16 t/(hm?.a), which is approximately 1.5 to 3.0 times of the global 
average of erosion (Mosaffaie and Talebi, 2014). Land use change makes the area susceptible to 
wind erosion, which has a significant effect on carbon dioxide emitted from the soil to the 
atmosphere due to affecting microorganisms' activity, soil moisture and temperature (Rahimi et al., 
2013). Wind erosion also causes some changes in soil moisture and temperature (Karmakar et al., 
2016). Soil can adapt to new conditions of land surface and change the soil organic carbon by 
mineralization in response to moisture and temperature (Chappell et al., 2016; Wang et al., 2017). 
In addition, erosion of surface soil by wind could have some effects on soil microenvironment, 
decomposition rate and root depth (Chappell et al., 2016). As previously mentioned, soil organic 
carbon storage is very sensitive to management practices and land use (Tan and Lal, 2005). Soil 
manipulation due to anthropogenic activities such as plowing, tillage and cultivation (Schlesinger, 
1984; Leiber-Sauheitl et al., 2013) increases the activity of microorganisms, erosion and high 
temperature because of mechanical activity and reduces moisture due to more evaporation 
(Joneidi et al., 2015), thus having an impact on carbon emissions rate. Therefore, recognizing 
these gas exchanges is essential for finding the right way of land management for more carbon 
from the atmosphere to the soil and reducing carbon emission from the soil to the atmosphere 
(Schwining and Sala, 2004; Wohlfahrt et al., 2008; Zeeman et al., 2010; Sadeghipour et al., 2015). 

The second result of the current study is that the highest and the lowest carbon dioxide 
emissions from the soil occurred during the warmest months of the year (summer) and the coldest 
months (winter), respectively. This carbon dioxide emission was only positively related to air 
temperature and had no significant relationship with soil moisture content. It can be because of the 
low moisture fluctuation in different seasons due to being located in a dry area. One of the most 
important factors affecting the activity of microorganisms and soil respiration is temperature. 
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Increasing the temperature to a certain extent increases the activity of microorganisms. Humidity 
can affect soil respiration rate (as an indication of high biological activity), mainly due to its effect 
on soil temperature. In the study area, the moisture of the soil is not affected for long periods, even 
during the humid months of the year, so that the effect of moisture on soil respiration and carbon 
dioxide emission is not severe. 

Many researchers reported that the highest and the least amounts of emission occur in summer 
and winter, respectively, which has a positive correlation with the temperature and soil biological 
activities (Inubushi et al., 2003; Ussiri and Lal, 2009; Santra et al., 2013). As reported by Fortin et 
al. (1996), increase in carbon dioxide emission with an increase in temperature is attributed to the 
increase in soil respiration. Changes in temperature, precipitation and evaporation will cause a 
significant change in organic matter turnover and carbon dioxide dynamic (Karmakar et al., 2016). 
As the temperature goes up, soil microorganisms and root respiration rise, which causes more 
emission. 

In addition, results represented that there is no significant difference between the control and 
erosion class of low, in term of carbon loss; hence reducing the amount of erosion by restoration 
and reclamation projects, even if result not in the complete elimination of soil erosion, can reduce 
the amount of carbon emission in significant quantities. 


5 Conclusions 


Accelerated erosion at watershed scale and the short-term timeframe can be considered as a source 
of carbon dioxide. Wind erosion will lead to a significant rising in carbon dioxide emission from 
the landscape into the atmosphere. Besides the wind erosion, increasing air temperature will also 
lead to an intensification of carbon dioxide emission. In this study, a significant interaction was 
observed between wind erosion and high temperature, demonstrating the synergetic effect of wind 
erosion and climate warming. Increasing carbon emission due to the wind erosion and its 
consequences will exacerbate global warming, in turn, will increase carbon emission in a 
repeating two-way manner. 

There are many factors involved in erosion including petrology, land geomorphology and 
topography, wind intensity, soil texture, vegetation cover, etc. That should be considered along 
with the biological properties of the soil. Finally, effect on soil respiration and carbon dioxide 
emission must be investigated, to be considered in management plans. It seems that due to the 
extensive soil degradation in Iran, especially in arid and semi-arid areas, the possibility of natural 
restoration is very low, and it is imperative to use mechanical reclamation techniques along with 
afforestation to prevent soil erosion. Thus, erosion control and reclamation of eroded soils can be 
important solutions for increasing soil quality and reducing the risks of greenhouse gases 
increasing. 

The positive relationship between carbon dioxide emission and air temperature indicates the 
high potential of global warming to increase carbon emission at a great scale. Global warming 
itself can lead to processes that increase soil erosion. Soil erosion also increases carbon emissions, 
and increases in GHG concentration, including carbon dioxide, causing global warming. 
Therefore, we are faced with a complex matrix of multilateral impact factors, which reveal the 
need to apply strategies to adapt to global warming. 
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